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We study the fluctuations of the directed polymer in 1+1 dimensions in a Gaussian random envi- 
ronment with a finite correlation length £ and at finite temperature. We address the correspondence 
between the geometrical transverse fluctuations of the directed polymer, described by its rough- 
ness, and the fluctuations of its free-energy, characterized by its two-point correlator. Analytical 
arguments are provided in favor of a generic scaling law between those quantities, at finite time, 
non-vanishing £ and explicit temperature dependence. Numerical results are in good agreement 
both for simulations on the discrete directed polymer and on a continuous directed polymer (with 
short-range correlated disorder). Applications to recent experiments on liquid crystals are discussed. 
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I. INTRODUCTION 

Brownian particles have provided in Physics one of the 
first example of systems whose geometrical properties dif- 
fer radically from those encountered in regular classical 
mechanics: instead of following 'smooth' (differentiable) 
trajectories, those particles follow continuous but 'rough' 
(non-differentiable) paths due to the persistent stochas- 
tic thermal forces they withstand at finite temperature. 
This phenomenon is described in statistical mechanics by 
a random walk, whose geometrical self-similarity at large 
scale is depicted by scaling laws and scaling exponents. 
These scalings are known to be distinctive features of uni- 
versality classes - gathering different physical phenomena 
sharing a common representation. A natural question re- 
garding such random paths pertains to the influence of 
the environment: how does a path in a uniform medium 
differ, for instance, from a path in a medium with random 
inhomogeneities? And how do such differences manifest 
in scaling properties? 

Beyond particle trajectories, those paths also describe 
generic interfaces or random manifolds between distinct 
phases. Example systems include imbibition fronts [1, 2], 
wetting and spreading interfaces [3, 4], cracks propagat- 
ing in paper [5], avalanches of pinned interfaces [6, 7], 
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burning fronts [8, 9], interfaces in magnetic [10-12] or 
ferroelectric materials [13, 14] and generic growth phe- 
nomena [15, 16]. A wide selection of those systems, al- 
though ranging from microscopic to macroscopic scale 
and presenting a large variety of microphysics, have 
been described as generic disordered elastic systems (see 
Refs. [17-23] for reviews). 

One class of random walk has received close attention 
in the past decades, the one-dimensional Kardar-Parisi- 
Zhang (KPZ) [24] universality class (see Refs. [25, 26] 
for recent reviews) , since it is related to questions of ex- 
tremely varied nature [27] ranging from Burgers equation 
in hydrodynamics [28] , directed polymers in random me- 
dia [29, 30] to the parabolic Anderson model [31], eigen- 
values of random matrices [32-34], vicious walkers [35- 
38], dynamics of cold atoms [39] and transport in ID 
stochastic [40] or deterministic models [41]. 

Here we address the link between the scaling proper- 
ties of the geometry of the directed polymer (through its 
scale-dependent roughness) and the scaling of its free- 
energy fluctuations. We examine in particular the role of 
temperature for those scalings and its possible interplay 
with finite disorder correlation length, especially at finite 
time. 

More precisely, let us denote by t the direction of time 
of the growing directed polymer (DP), which also repre- 
sents the 'scale' at which a generic interface is examined, 
and by y the transverse direction in which fluctuations 
occur (see Fig. 1). One important aspect of the KPZ 
universality class is that the fluctuations of the (suitably 
centered) free-energy F(t,y) are expected to behave as 
follows in the large time asymptotics: 



F(t,y) 



(*- 



,i ( y 



(i) 



where a and b are constants which depend on the phys- 
ical parameters of the system (e.g. temperature T, elas- 
ticity c, disorder strength D and disorder correlation 
length £ in the model we use, see section II), and x(y) 
is a ^independent 'random variable' whose distribution 
characterizes the fluctuations of y-dependent observables. 
Note that in some systems, e.g. liquid crystals, the quan- 
tity F(t,y) scaling as (1) is, up to a translation, the 
height of an interface and not its free-energy [42-44] . For 
a directed polymer with fixed endpoints (as depicted in 
Figs. 1 and 2) and with uncorrelated random environ- 
ment (£ = 0) it has been shown [45, 46] that x(y), consid- 
ered as a (stochastic) function of y, is equivalent in distri- 
bution to the so-called Airy2 process (minus a parabola) , 
which is a stationary determinantal point process. It fol- 
lows from this relation (1) that, in the large t limit, the 
fluctuations of F(t,y), at fixed y, are described by the 
Tracy- Widom distribution Ti [47]. The latter distribu- 
tion J-~2, which can be written as a Fredholm determinant 
involving the Airy kernel, describes the fluctuations of 
the largest eigenvalue of random matrices belonging to 
the Gaussian Unitary Ensemble (GUE). Still at £ = 0, 
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FIG. 1. Schematic view of a continuous directed polymer of 
trajectory y(t) starting in y — at time and arriving in yi at 
time ti, in a quenched random potential V(t, y). These point- 
to-point configurations correspond to the "droplet geometry" 
of growth models. 



remarkably, it was shown that at finite time t, the dis- 
tribution of F(t, y), for fixed y, can still be written as a 
Fredholm determinant, which involves a non trivial finite 
t generalization of the Airy kernel [48-51]. 

In this paper we focus on the effects of finite temper- 
ature and finite disorder correlation length (£ > or 
discrete DP) at finite time for which we propose a gener- 
alized scaling relation which reads formally 



F(t,y) ~ak/B(i) 



V 



(2) 



where B(t) is the roughness of the directed polymer, de- 
fined as the second cumulant of the transverse fluctua- 
tions along y at fixed time t. A more accurate statement 
is made in section III B in terms of the two-point corre- 
lator of F(t,y), see Eq. (22). The large time asymp- 
totics (1) is then recovered from B(t) t ~?2° . Our 
motivation for investigating the case £ > comes from 
the study of physical or chemical experiments where the 
disordered potential always presents a finite correlation 
length, which, albeit microscopic and often inaccessible 
to direct measurement, may still induce observable effects 
at large scales [23]. In particular, we probe the depen- 
dence of the constants a, b on the system parameters by 
distinguishing a high- and a low-temperature regime due 
to £ > 0. The appearance of those regimes generalizes 
the ones affecting some specific observables (such as the 
roughness of the interface, studied in Refs. [52, 53]). 

The directed polymer model we study is defined in sec- 
tion II and the generalized scaling form we propose is 
presented and discussed analytically in section III. Nu- 
merical results are gathered in section IV for the discrete 
DP and section V for its continuous version. We discuss 
experimental implications and draw our conclusions in 
section VI. The details of some computations are gath- 
ered in the appendices. 
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FIG. 2. (Color online) Geometry of the discrete DP model 
we consider. Allowed steps are y(t + 1) — y(t) — ±1. The 
longitudinal and transverse axes correspond to the directions 
t and y respectively. 



II. MODEL AND SCOPE OF THE STUDY 

A. DP formulation 

We focus on the directed polymer formulation: the tra- 
jectory of the polymer is described by a continuous coor- 
dinate y(t), starting from y — at £ = (see Fig. 1), and 
growing in a random potential V(t, y) along direction £. 
The energy of a trajectory y(t) of duration t\ is the sum 
of elastic and disorder contributions: 

H[y, V; t x ] = J 1 dt { \ [d t y(t)j 2 + V (t, y(t)) } . (3) 

The first term flattens the polymer by penalizing its de- 
formations (with an intensity encoded in the elastic con- 
stant c), while the second term tends to deform it. At 
fixed disorder V, the weight of all trajectories starting 
from and arriving in yi at time £i is given by the path 
integral 

Zv(h,yi) = f yitl) ~ Vl Vy(t) e-^lvy-M] ) (4) 

Jy(Q)=Q 

(we set the Boltzmann constant equal to 1 and denote 
thereafter the inverse temperature by /3 = yp ) . The mean 
value of an observable O depending on the value of y at 
final time £ reads, at fixed disorder 



(0[y(t)]h 



JdyO(y)Z v (t,y) 
JdyZ v (t,y) 



(5) 



Here and thereafter, the integrals J dy over the DP end- 
point run by convention on the real line. We also consider 
a discrete version of the same system, illustrated in Fig. 2 
and described in section IV. 

The weight Zy(t,y) is not normalized to unity 
(/ dy Zv(t,y) ^ 1 in general) but the path integral (4) 
is chosen so that at zero disorder J dy Zv=o(t,y) = 1 at 
all times. With this choice of normalization, it is known 
from the Feynman-Kac formula [54-57] (see also Ref. [58] 



for a discussion) that the weight Zy (£, y) evolves accord- 
ing to the 'stochastic heat equation' [59-61] 

d t Z v (t,y) = ^d 2 v Z v (t,y) - ±V(t,y)Z v (t,y) , (6) 



2c 



T 



while the free-energy Fy{t,y) = — T\og Zy(t, y) obeys 
the Kardar-Parisi-Zhang equation [24, 56] 



d t F v (t,y) 



2c 



[d y F v (t,y)] 2 + V(t,y) 



(7) 

The presence of the so-called non-linear 'KPZ term' 
makes the study of this equation particularly difficult; in 
particular the typical extension of the excursions along 
the direction y does not scale diffusively at large times 
(i.e. y ~ £2) but superdiffusively (y ~ £*» where £ = | is 
the KPZ exponent) [24, 28, 32, 45, 56, 62]. 

The distribution of the random potential V(t,y) de- 
termines the statistical properties of the free-energy. We 
denote by an overline (£%(£)] )y the statistical average 
over the disorder V. One case is well understood: when 
V(t, y) is a centered Gaussian uncorrelated random po- 
tential (£ = 0), that is, when V(t,y) has a Gaussian 
distribution with zero mean and two-point correlator 



V(t,y)V(t',y')=DS(t'-t)S(y'-y) 



(8) 



it is known for long [56] that the infinite-time distribution 
of the free-energy is that of a two-sided Brownian walk: 



Prob[F(y)] oc exp 



1 T 
2cD 



dy[d v F(y)\ 



(9) 



In other words, Fy(t, y) has a steady-state Gaussian dis- 
tribution whose correlator reads: 



[F v (t,y') - F v (t,y)Y 



cD 
~T 



W y\ 



(10) 



The value | for the KPZ exponent is known since 
the work of Henley, Huse and Fisher [56] and Kardar, 
Parisi and Zhang [24] and Kardar [62, 63] but this re- 
sult has been proven in mathematical framework much 
later [32, 45]. The free-energy distribution has recently 
been determined by a variety of methods both in the 
physics [48, 51, 64] and in the mathematics [49, 50, 65] 
communities (see Ref. [66] for a review). 



B. Generalization of the free-energy fluctuation 
scaling 

In the large but finite time limit, Prahofer and 
Spohn [45] have shown that at £ = the correlator 



C(t,y) = [F v (t,y)-F v (t,0)Y 
obeys the following scaling relation 

C(t,y) = ti y PS (yt~i) , 



(11) 
(12) 
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with gps(y) ~ %jr\y\ for small |y| (which thus gives (10) 
in the limit t — > oo) while gps(j/) saturates to a con- 
stant at large \y\. This function gps(y) is> u p to non- 
universal longitudinal and transverse length scales, the 
mean square displacement of the Airy2 process discussed 
below in Eq. (42). 

Our aim in this paper is to propose and test an exten- 
sion of the scaling relation (12) (i) at finite time and/or 
(ii) for a disorder V(t,y) presenting short-range correla- 
tions. Much less is known in those two cases. To imple- 
ment short-range correlations (on a transverse scale £) in 
the disorder distribution, we assume that it is still zero- 
mean Gaussian distributed with correlations of the form 



V(t, y)V(f, y>) = D5{t' ~ t)R £ (y' - y) 



(13) 



Such correlations in direction y on a range of order £ are 
described by a rounded delta peak R.^(y), normalized to 
unity (f dy R^(y) — 1) and by the strength D of the dis- 
order. Note that even for £ > the distribution of V is 
invariant by translation along y, as in the uncorrelated 
case (8) (which corresponds to £ = 0). For simplicity we 
may assume in explicit examples that (y) is a normal- 
ized Gaussian of variance 2£ 2 : 



R 



Gauss 



(y) = 



i 



VW 2 



(14) 



Note that the knowledge of the fluctuations of the free- 
energy Fy is not sufficient to retrieve the mean value of 
physical observables from (5): in general the full distri- 
bution of Fy is required. The free-energy fluctuations 
described by the correlator C(i, y) still provide physically 
relevant information, for instance combined to scaling 
arguments [67], or as the starting point of 'toy-model' 
approach where the free-energy distribution is approxi- 
mated to be Gaussian [53, 68-71] (see also part IIIC). 



III. SCALING OF THE FREE-ENERGY 
CORRELATOR C s (t,y) 
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FIG. 3. Schematic plot of the correlator C$(t,y) of F, as a 
function of y at fixed t. It is rounded close to the origin (\y\ < 
£), behaves in D\y\ for intermediate length (£ < \y\ < it) and 
goes to a constant at large length (\y\ > It). 



where c|r 



t j Q g 2^rt j g ^ e e i as tj c contribution, which 
fully captures the initial condition, while Fy(t,y) is in- 
variant by translation along y in distribution (in math- 
ematical terms [49] F v (t,y) is 'stationary' along y). It 
represents the contribution of the disordered potential to 
the free-energy (Fy(t,y)\y=o = 0). Note also from (15) 
that the initial condition Zy(0,y) = S(y) simply writes 
Fy(0, y) = 0. This form of initial condition is technically 
different from the 'sharp wedge' often considered to pin 
the polymer in y = at initial time [50] and which also 
ensures Zy(0,y) = S(y). 

This translational invariance allows to define our cor- 
relator of interest: 



C i (t,y'-y)= [F v (t,y>) - F v (t,y)] (16) 

where we have made explicit the dependence in £ arising 
from the distribution of V. The STS ensures that (16) 
depends in y and y' only through the difference y' — y, 
while the invariance of the distribution of the disorder 
V by the reflection V H> — V ensures that the function 
C^(t, y) is even with respect to its argument y. Note that 
at £ = 0, the steady-state result (10) implies 



A. A generalized scaling relation 

Let us first identify a definition of the frcc-cncrgy cor- 
relator which is suitable to study its finite-time scal- 
ing. The correlator [Fy (t, 7/2) — Fv (t, 2/i)] 2 is invariant 
by translation along y only in the infinite time limit (10). 
To extend this property at finite time, one may take ad- 
vantage of the 'statistical tilt symmetry' (STS) [72-74] of 
the model which ensures that the free-energy splits into 
two contributions (see e.g. Ref. [58] for a derivation at 
non-zero £): 



, , V T, 2nTt - . . 

Fv(t,v) = c- + -log + F v (t,y 15 

It 2 ^ c 

clastic contribution disorder 

contribution 



cD 



(17) 



The correlator C^(t,y) is actually the connected correla- 
tor of the full free-energy Fy(t,y) (see appendix A) 



C i (t,y'-y)=[Fy(t,y')-Fy(t ) y)}' 



F v (t,y')-F v (t,y) 



(18) 



The expression (16) is useful to explicit the translational 
invariance while (18) enables to consider cases where no 
decomposition such as (15) is available (as for the dis- 
crete DP). 

How can we interpret the crossover from the ini- 
tial condition C\(0, y) — to the steady-state result 
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Q(oo,y) oc \y\1 The weight cx e ~P^~fSF v (t, y ) of a tra _ 
jectory ending in y at time t depicts a particle of position 
y (time t being fixed) in an "effective potential" made up 
of a parabolic potential describing thermal fluctuations 
and of an "effective disorder" F, which recapitulates the 
disorder landscape V perceived by the polymer from its 
starting point. Heuristically, one expects that F remains 
almost flat at initial times since the polymer end-point 
y(t) has not explored much of its random environment, 
while at larger times F becomes a "random force poten- 
tial" (with cuspy correlator C^(t,y) cx |y|), at least on a 
transverse region \y\ < £ t of typical size it in which the 
polymer endpoint has mainly been confined. 

We show in appendix B that this intuitive picture is in- 
deed correct: at all finite times, lim^i^^ d y C^(t, y) = 0. 
The correlator C^(t, y) thus has to switch from the abso- 
lute value behavior \y\ for |y| < £ t to a plateau for \y\ > £ t 
(see Fig. 3), at some scale i t increasing and diverging with 
t. We assume that the effect of the correlation length £ is 
to round C^(t,y) at small \y\ < £. The following scaling 
relation is thus expected to hold: 



C e (t,y) = ttCt-^y) 



(19) 



where the scaling function C^{y) depends on the physical 
parameters c, D, T and £. At zero £, it is compatible with 
the large time behavior (12), provided that it ~ for 
t — > oo. We now have to identify the crossover length i t 
at non-zero £ and finite time. 



B. Asymptotic transverse fluctuations 

The variance of the DP endpoint y(t), called the rough- 
ness, is the simplest length quantifying the spatial exten- 
sion of the polymer at a given time t: 



B(t) = {y{ty)v 



Jdyy 2 Z v (t,y) 
JdyZ v (t,y) 



(20) 



It is known that the roughness presents at small times a 
diffusive regime (£ t h = \) and at large times a random 
manifold (RM) superdiffusive regime (Crm = §) [24, 32, 
56]. In terms of power laws of the time t, we have the 
two asymptotic regimes 



(21) 



(see also Refs. [53, 58] when the disorder correlation 
length £ is nonzero). We propose that i t ~ \J B(t) 
in (19), in other words 




(22) 



Before testing numerically this scaling law in different 
models (see sections IV and V) , we discuss analytical ar- 
guments in favor of (22). On one hand, linearizing the 




FIG. 4. (Color online) Graph of the correlator Cf n {t, y) of F 
obtained in the linear approximation (C23), as a function of y 
for different times t (parameters are c=l, D = l, T=l and 
£ = 0.3). Time increases geometrically from bottom (t = 2~ 4 ) 
to top (t = 2 10 ) curves. 



dynamics at short time in Eq. (7) one finds that the fluc- 
tuations are diffusive: B(t) = — . In this approximation, 
the correlations of the free-energy rescale as follows (see 
Appendix C): 



Cj n {t,y) = 



cD 
~T 



y/Bth(t) C h 



(23) 



Ti 



is the thermal roughness. The scaling 



where B t h(t) 

function C^ n (y), given in (C27), function of the properly 

adimensional variables y and £, is independent of the 
parameters c, D and T. The behavior of C^ n (t,y) is 
as qualitatively expected with a rounding on a scale £, 
a plateau at large y and a developing linear behavior in 
between, as plotted in Fig. 4. 

On the other hand, at zero £ and large time, the re- 
sult (12) of Prahofer and Spohn is compatible with (22): 
this corresponds to the RM asymptotics B(t) ~ ts . To 
address the RM regime at non-zero £, we study in the 
next paragraph a toy-model approach of the DP. 



C. A winged DP toy model 

In the infinite time limit and at zero £, the free-energy 
Fv(t,y) is translationally invariant in distribution and 
has a Gaussian distribution (9). The statistical tilt sym- 
metry expresses that the free-energy splits at finite time 
in two contributions (15), of which only Fy(t, y) is trans- 
lationally invariant. The idea of the toy model [68-70, 74] 
is to assume that for large t and finite £ the distribution 
of the reduced free-energy (denoted F toy (t,y)) remains 
in a good approximation Gaussian, with zero mean and 
correlations 



Ct° y (t, y'-y)= [F*°y(t, y>) - F*°y(t, y)] : 



(24) 
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It has been shown [53] that the roughness B(t) can be 
computed in a Gaussian Variational Method (GVM) ap- 
proximation, for a correlator C^ oy (t,y) of the Fourier 
form 

Ct° y (t,y) = J g^[l-cos(A 2/ )]i?i oy (t,A) , (25) 

with i?~ oy (i, A) = De~ x2 t 2 . For f = 0, this writing yields 

the £ = infinite time-result (10) C(t,y) = D\y\, pro- 
vided that D = For f > 0, i?|° y (t, A) encodes for F 
the finite correlation length of the disorder V: the ab- 
solute value becomes rounded around the origin up to a 
scale ~ £. The problem of the form (25) of the correlator 
is that it implies C^ oy (t,y) ~ D\y\ at large \y\, while we 
have seen that C^(t,y) goes to a constant in the limit 
\y\ — > oo at all finite times t. 

To overcome the discrepancy between this exact re- 
sult and the model i?| oy (i, A) = De~ x t , one may rather 

study a correlator C^ oy (t,y) of the form (25) with addi- 
tional saturation 'wings' 

Rfy { t,X)=D^^e-^ , (26) 

which still presents a behavior C toy (t,y) ~ D\y\ at in- 
termediate y (£ -C | y | -C 4)) but goes to a constant at 
large \y\ (\y\ 3> ^t). Note that more generically one can 
consider a toy correlator of the form 

Rf y (t 7 \)=Df 1 (\£ t )f 2 (\0 , (27) 

in which the rounding due to the finite disorder correla- 
tion length appears in the factor /2(A£) and the crossover 
at scale £ t in the function /i (X£t) [75] . In particular, scal- 
ing properties of the roughness arise from the form (27). 
Note that this form makes sense in principle only when 
the scales £ t and £ are separated enough (with £ t > £). 
We compute in appendix D the roughness of this model 
in the Gaussian Variational Method (GVM) approxima- 
tion, in the large £ t limit. The computation yields, in the 
random manifold regime 

B(t) = B(t) |,_ 1=0 - ^-t 2 + 0{£~ t 2 ) , (28) 

where the RM roughness in the absence of wings (£^=0) 
is given by 

For the scale £t not to destroy the | exponent of the RM 
regime, one must have £ t growing at least as 

4 ~(f)M OO) 

In other words, the wings of C(t, y) have to appear at a 
scale larger than \/ B{t) in the RM regime. 



IV. NUMERICAL SIMULATIONS AT HIGH 
TEMPERATURE: THE DISCRETE DP 

A. Model 

We use numerical simulations of the DP model in order 
to test the scaling properties of free-energy fluctuations 
at finite temperature. Following the geometry described 
in Fig. 2 we perform numerical simulations of a discrete 
DP model [63] with the solid-on-solid (SOS) restriction 
\y(t + 1) — y(t)\ = 1. A site-dependent zero-mean uncor- 
related Gaussian disorder potential V t . y of intensity D is 
used: 

V t ,yV V ,y. = D6t, V &y,y> 

The energy of a given configuration of the DP is given by 
the sum of the site energies along the path y(t). 

Given a disorder realization characterized by Vt, y , the 
weight Z ttV of a polymer starting in (0, 0) and ending in 
(t, y) is given by the following recursion: 

Z t , y = e~ PVt - v (Zt-i, v -i + Zt-i,v+i) (32) 

with zero time initial condition Z y = 5o >y . Therefore, 
the probability to observe a polymer ending in (t, y) is 
Zt,y/^2 y , Zt,y'- Due to the recursion relation, Z tlV grows 
exponentially with time t. To avoid numerical instability, 
all weights Z t y at fixed t are divided by the largest one, 
which does not change the polymer ending probability. In 
terms of the weight Zt, y the free-energy of the polymer 
starting in (0, 0) and ending in (t, y) is therefore defined 
as 

F '-- Th {A)' (33) 

and the free-energy fluctuations are measured in terms of 
the connected correlation 

C(t,y) = [F t , v - F^ y ,] 2 - [F t , y - F w f , (34) 

As discussed in paragraph III A, this definition is equiv- 
alent in the continuum to that involving F(t,y) in (16). 
The index £ is dropped from C(t, y) since the disorder 
is uncorrelated as in (8), and since no lengthscale be- 
low the lattice spacing can be considered in this discrete 
DP model. The roughness is defined as the mean square 
displacement of the free end, as in (20) for the continu- 
ous DP 

y i-V 

Finally, note the correspondence between the discrete 
and continuous parameters [see (E8) with lattice spacings 
a and b taken to 1)]: the continuum model (6) which is 
the limit of the discrete model (32) in the large size limit 
has the parameters 

temperature T, elasticity c = T, disorder D (36) 
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FIG. 5. (Color online) Bare free-energy fluctuations C(t,y) 
for the discrete DP at finite temperature T — 4 and for two 
different times t. The late increase of C(t,y) for t = 512 is a 
finite size effect related to the solid-on-solid restriction of the 
DP model. 
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FIG. 7. (Color online) Rescaled free-energy fluctuations 
C(t, y) for the discrete DP at large times and high tem- 
peratures, according to the scaling (40). The inset 
shows the bare data C(t, y) against y for timescales t = 
4096, 8192, 16384, 32768 from bottom to top. Thin (thick) 
lines correspond to T = 4 (T = 8). 




FIG. 6. (Color online) Rescaled free-energy fluctuations 
C(t,y) for the discrete DP at finite low temperature T = 0.5 
and various times, according to the scaling (37) with £ = 
Crm = §- The inset shows the bare data C(t, y) against y for 
timescales t = 4096, 8192, 16384, 32768 from bottom to top. 

thus providing a comparison of the temperature depen- 
dence of the two models through the correspondence 
c = T. A generic study of the continuum limit of the 
discrete DP is presented in Appendix E. 



B. Results 

In what follows, results are presented for the discrete 
DP model with D = 1 and disorder averages performed 
on 10 disorder realizations. Fig. 5 displays typical curves 
for the free-energy fluctuations C(t, y) at a temperature 
T = 4. Two different time scales are presented, t = 512 
and t = 8192, in order to show how the finite size of 
the system affects free-energy fluctuations. It is clear 
that one observes for both timescales the characteristic 



\y\ behavior at small transverse lengthscales and then a 
crossover to saturation at larger timescales. In addition, 
when the polymer length is small an increase of C (t, y) is 
observed for large y, as can be observed for the t = 512 
data. This last point is a finite-size effect related to the 
fact that, in contrast to the continuous case, for the dis- 
crete DP model used here the polymer endpoint is con- 
strained at all times to remain in the cone \y{t)\ < t, see 
Fig. 2. This finite-size effect only plays a role when ana- 
lyzing small DP lengths and in the following we discard 
this finite-size regime in order to better compare different 
curves. 

We now test the different scaling properties for C (t, y) . 
With respect to time t and in the large time limit [45], 
C(t, y) is expected to scale according to 

C(t,y)~t<d(£), (37) 

with £ = £ R m = | at large t, which corresponds to 
the scaling relation in Eq. (12). This behavior has 
been reported in Ref. [69] for a finite low temperature 
value T « 0.14 (j3 — 7) and moderated polymer sizes 
t = 512, 1024. In Fig. 6 we test the scaling given in 
Eq. (37) for temperature T = 0.5, which is larger than 
the one used in Ref. [69] but still in a low temperature 
regime (see below). As shown in the figure the scaling 
works satisfactorily for the large time scales used. This 
result confirms the scaling probed in Ref. [69] and subse- 
quently analytically obtained in Ref. [45]. 

In order to go beyond the large time scaling of Eq. (37) , 
we now discuss the temperature dependence of free- 
energy fluctuations with respect to the roughness. It is 
first important to properly consider the finite tempera- 
ture scaling of the roughness, which has been discussed 
in Refs. [52, 53, 76]. In the high temperature regime 
and within the continuum limit of the DP, the roughness 
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FIG. 8. (Color online) Rescaled free-energy fluctuations 
C(t, y) for the discrete DP at large times and low tem- 
peratures, according to the scaling (40). The inset 
shows the bare data C(t, y) against y for timescales t = 
4096,8192,16384,32768 from bottom to top. Thin (thick) 
lines correspond to T = 0.5 (T = 1). 



scales asymptotically as 



B(t) 



(t « L r ) 



ith Lt = 



(38) 

where Cth — 1/2 is the thermal roughness exponent 
and Crm = §■ The temperature dependence of the 
roughness is described at large scale by the thorn expo- 
nent p = — 1/3. The thermal lengthscale Lt separates 
short-scale thermal fluctuations characterized by Cth 
from large-scale disorder-induced fluctuations of expo- 
nent Crm- It has been shown [76] that in the high 
temperature regime the thermal length scale grows as 
Lt ~ T 1 /° F /(cD 2 ), where F = 1/5 is the Flory expo- 
nent [52, 53]. In the discrete version of the DP one has 
c = T [see (36) and Appendix E] and therefore the scaling 
of the roughness (38) now reads 



B(t) 



t 2(th 

T 2 \ 2 i> 



(t « L T ) T 4 
(£)^hm {t>>LT) withL T -- 



(39) 

In this last case, remarkably, the short-time behavior is 
temperature-independent while the large time prefactor 
scales with temperature as T 41? = 7 1 - 4 / 3 . 

These scaling properties of the roughness can be di- 
rectly incorporated into the scaling properties of free- 
energy correlations using the whole roughness to define 
the characteristic transverse scale as in the scaling rela- 
tion (22) which, in the absence of £, writes 



(40) 



First note that the inset of Fig. 9 shows how the sat- 
uration regime of free-energy fluctuations is reached at 
smaller transverse lengthscale y when increasing the tem- 
perature. If, as suggested by the scaling relation (40), the 
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FIG. 9. (Color online) Rescaled free-energy fluctuations 
C(t, y) for the discrete DP at a fixed large time t = 16384 
and from low to high temperatures, according to the scal- 
ing (40). The inset shows the bare data C(t,y) against y. 
Temperatures are T = 0.5, 1, 2, 4, 8 from top to bottom. 




FIG. 10. (Color online) Full temperature-independent 
rescaled free-energy fluctuations C(t,y) for the discrete DP 
at a fixed large time, according to the scaling (41). Same 
data as in Fig. 9. The inset shows the parameter D of (41) 
as a function of T. 



crossover is dictated by the time dependent roughness 
B(t), this is in agreement with the negative value of the 
thorn exponent ruling the temperature dependence of the 
roughness [52, 53, 76]. The main panel of Fig. 7 illustrates 
the validity of the scaling relation (40) that incorporates 
the temperature dependence through the roughness func- 
tion, in the range of high-temperature roughness. 

Since the temperature dependence is included in the 
roughness B(t), this scaling relation can also be probed 
at lower temperatures. This is done in Fig. 8, were data 
corresponding to two low temperatures T = 0.5, 1 and 
different time scales are collapsed on a single curve. How- 
ever, the high- and low-temperature regimes are not nec- 
essarily described by the same rescaling function C(y). 
Fig. 9 actually shows, for a single value t — 16384, that 
although high- and low-temperature data collapse on a 
single curve proper to each regime, there is a crossover 
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FIG. 11. (Color online) Comparison between C 1 (y) denned in 
Eq. (41), extracted from the numerical data shown in Fig. 10 
and Agps(y/B) (dashed line), where <?ps(y) is denned in 
Eq. (42) while A and B are two fitting parameters (deter- 
mined to A ~ 4.08 ± 0.05 and B ~ 2.07 ± 0.06). The numeri- 
cal evaluation of the exact expression for pps (y) was done in 
Refs. [77, 78]. 

between these two limiting cases. 

As suggested by the scaling relation (23), the slope 
of the linear initial growth of the free-energy correlator 
should fully account for the temperature dependence. As 
extensively discussed in [58] , it can in fact be argued that 
the (D, Independence of the scaling law (22) can be ab- 
sorbed in a single pref actor D, which in the case of the 
discrete DP model reads 

C(t,y)=D^W)C 1 (-^), (41) 

where D — D(D,T) and the function C' 1 (y) is indepen- 
dent of the parameters D and T. Indeed, the prefactor D 
in (41) is the slope of the correlator C(t, y) in the regime 
\y\ ~ \J -^W where C(t,y) cx \y\. It is known from the 
large-time limit at zero £ that D = see (17) (and 
also (C25) in the diffusive regime). In the discrete DP 
case (c = T), D is thus expected to be temperature in- 
dependent in the high temperature regime, D = D, and 
to behave as T 2 / 3 in the low temperature regime [58]. 
In order to test this, we rescale the data in Fig. 9 onto a 
single universal master curve as shown in Fig. 10. To this 
purpose, we fix D = D = 1 for T = 4, 8 and ensure the 
collapse of low temperature correlators C (t, y) by propos- 
ing values for the parameter D/D for T = 0.5, 1, 2, which 
are plotted in the inset of Fig. 10. Although the value 
of D is decreasing with T we did not observe the T 2 / 3 
behavior, possibly because we are not reaching the corre- 
sponding low temperature asymptotic regime for the dis- 
crete DP model. However it is remarkable to observe the 
crossover towards the lower-temperature regime, which 
reflects the influence of the lattice spacing even though 
no lengthscale is defined below it. 

In summary, we have shown in this section that numer- 
ical simulations of the discrete version of the DP model 
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FIG. 12. (Color online) Graph for the continuous DP of the 
rescaled correlator C TCBC (t,y) defined in (45) as a function 
of y = y/WB(t) for different times t, at low temperature. 
Parameters are c = 1, D = 8, T = 0.4 and £ = 2. Disor- 
der average is performed over 1092 realizations. Inset (a): 
Roughness B(t) in log-log scale with the thermal (green line) 
and random manifold (red line) asymptotic regimes. Inset 
(b): Corresponding bare correlator C(t,y). Time varies from 
35 to 135 by steps of 25 from bottom to top curves. 



with the SOS constraint are fully compatible with the 
scaling arguments suggesting that \J B{t) is the relevant 
transverse scale for free-energy fluctuations. Moreover, 
we have also shown that all the temperature dependence 
can only be accounted for by using, in addition to the 
proper transverse length scale yj B(t) the scale of linear 
free-energy fluctuations at small y, given by the func- 
tion D. This universal behavior being established, one 
would naturally expect from the result of Refs. [45, 46] 
in Eq. (12), that this correlator C' 1 (j/) in Eq. (41) can be 
expressed in terms of the mean square displacement of 
the Airy 2 process, ^4.2 (z/) ? namely 



C'iy) = Ag PS (y/B) , g PS (y) = ([A 2 {y) - .4 2 (0)] 2 ) , 

(42) 

where A and B are longitudinal and transverse 
parameter-dependent lengthscales. In Fig. 11 we test 
this relation (42) by adjusting the parameters A and B 
using least square fittings and we find indeed a very good 
collapse of our numerical data and the exact expression 
°f ffps(y) obtained in Ref. [45]. Note that the numeri- 
cal evaluation of <?ps(y) was done in Ref. [77, 78], where 
precise numerical techniques were developed to compute 
Fredholm determinants with high precision. 

In the following section we extend this analysis to lower 
temperatures and short-range correlations by studying a 
continuous version of the DP model. 
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FIG. 13. (Color online) Same graph as Fig. 12 but at higher 
temperature T = 1.5 (disorder average is over 446 realiza- 
tions) . 



V. NUMERICAL SIMULATIONS AT LOWER 
TEMPERATURE: THE CONTINUOUS DP 

A. Numerical approach 

To probe the implications for the scaling laws of short- 
range correlations (£ > 0) of the disordered potential 
V(t,y), we have simulated a continuous version of the 
directed polymer. Directly integrating the partial differ- 
ential equation (PDE) (6) for Zv{t,y) is difficult since, 
especially at low temperature, the weight concentrates 
exponentially in the most favorable regions for the poly- 
mer, preventing the whole space to be embraced. Taking 
the logarithm and considering the PDE (7) for Fv(t, y) is 
also problematic because of the singular initial condition 
Z v (0, y) = 8(y). We took advantage of the STS (15) by 
directly simulating the PDE for the reduced free-energy 

Fv(t,y) 

d t F v (t,y) = ±-%F v (t,y) - ^[d y F v (t,y)] 2 

- V -d v F v (t,y) + V{t,y) (43) 

whose initial condition is simply Fy(0,y) — 0. Another 
benefit of considering F v (t,y) instead of F v (t,y) is the 

2 

removal of the quadratic contribution c|r which eclipses 
the disorder contribution in (6) at short time. 

The correlated disordered potential V(t, y) is con- 
structed as follows: independent random variables Vij 
are drawn from a centered Gaussian distribution of vari- 
ance D gIld on a grid of coordinates (t, y) — (j£t,j£y)- 
The continuum V(t, y) is defined as the two-dimensional 
cubic spline of the on the grid. Its distribution is 

Gaussian and characterized by its two-point correlator, 
which can be analytically computed [58] and takes the 
form 



V(t,y)V(t',y') = £>%(*' - t)R iy (y' - y) , (44) 



FIG. 14. (Color online) Plot of the rescaled correlator 
D~ 1 C' rcsc (t, y) as a function of y at fixed large time t and 
for different temperatures T £ {.4, .55, 1.2, 1.5, 1.8}. Inset 
(a): corresponding original correlator C(t,y). Lighter color 
corresponds to lower temperature. Inset (b): D as a func- 
tion of T (up to a numerical factor) in log-log scale. Points 
are determined by least-square minimization (see text) with 
respect to a reference curve at T — 2 which fixes D|t=2 = 1. 
The line gives the expected slope of the large temperature 
asymptotics D oc h . 

where is a normalized smooth delta function of 
width £. The amplitude of disorder is D = D slld t;t£,y 
and our correlation length of interest in the transverse 
direction is £ = £ y . In numerical simulations, we took 
& = 1, £y = 2, L> grid = 4 (and thus £ = 2 and D = 8). 
Simulations are run on a finite window [— y m ,y m ] in the 
transverse direction, outside of which the disorder po- 
tential is set to 0. Besides, since the equation (43) is ill- 
defined at t = 0, simulations were run starting at small 
initial time to with thermal initial conditions. A complete 
presentation of the numerical procedure can be found in 
Ref. [58]. 

B. Results 

For numerical simplicity, the scaling (22) was tested 
graphically at fixed £, defining a rescaled transverse co- 
ordinate y = y/^jB{t) and a rescaled correlator 

C rcsc (t, y) = -jL=Ct (t, yVW)) (45) 

Indeed testing the full scaling (22) where £ is also rescaled 
by \J B(t) would imply to measure C^(t, y) for many dif- 
ferent values of £ which was not numerically accessible 
[see below for an analytical explanation on the rescaling 
\/y/B(€j of (22)]. Note that in the discrete DP it was 
sufficient to study the scaling function of (40) without 
taking £ into account (no length below the lattice spac- 
ing, which plays the role of £, can be considered). 

Low-temperature results are shown on Fig. f2: the 
curves of (t, y) at different times superimpose upon the 
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rescaling (45) , with a slight discrepancy around the origin 
due to the fixed disorder correlation length £. At higher 
temperature (Fig. 13) this scaling remains valid, the ef- 
fect of the finite £ being as expected less important. In 
the evaluation of C Teac (t, y) from (45), the roughness B(t) 
is determined from the numerical results for F(t, y) us- 
ing (5) to evaluate the thermal average (y(t) 2 )v a t fixed 
disorder, averaging afterwards over the realizations of dis- 
order. For completeness, the graphs of the roughness are 
given in the inset (a) of Figs. 12 and 13. 

Similarly to the case of the discrete polymer (41), it can 
be argued that the (c, D, Independence of the scaling 
law (22) can be absorbed into a single prefactor D [58] 



C ( (t,y) = D^W) C x _ 



(46) 



where the function C^(y) does not depend on (c,D,T). 

By tuning appropriately D, we show on Fig. 14 from the 
numerical evaluation of the correlator that this is indeed 
the case. The function D(T) at fixed c and £ is evalu- 
ated by fixing a reference curve C vcsc (t,y) at T = T rc f 
(Tref = 2 in Fig. 14) and finding the best D(T) which 
minimizes the distance between C Tesc (t,y) and the refer- 
ence curve using the least square method. We checked 
that the result does not depend on the choice of T rc f nor 
on the choice of the fixed t, within numerical uncertainty. 
This method allows a determination of D(T) up to a T- 
independent constant which imposes D(T Te {) = 1. We re- 
fer the reader to Ref. [58] for an in-depth analytical and 
numerical study of D, regarding its scaling and physical 
interpretation. As displayed in the inset (b) of Fig. 14, 
the measured D is compatible with the predicted high T 
behavior D = ^? and a saturation at low T. 

We can actually provide a short argument e xplain ing 
the rescaling of the correlation length £ by yj B(t) ap- 
pearing in (46) , that we could not probe numerically. As 
the rounding of C(t, y) is quadratic in y for small y (see 
also (C21)), we expect that Cl(y) scales as 

&\{y) a bt for \y\ « . (47) 



The rescaling f = £/y/B(t) in (46) ensures thus that, in 
the large time limit of (46), C^{t, y) remains independent 
of t at small y: 



=, , . 3/rsO D 2 
t->oo 2t;y/ir 



(48) 



Without the precise rescaling £/y/B(t) in (46) this ex- 
pression would become singular. 

Note that because of the rounding at small y the fit- 
ting procedure of C' 1 (i, y) with respect to gps, as done 
for the discrete DP in Fig. 11, cannot be performed since 
the Airy 2 process is not adapted to a finite disorder cor- 
relation length £. 



VI. DISCUSSION AND CONCLUSION 

In the one-dimensional KPZ universality class, a suit- 
ably centered observable F(t, y) presents, in the large 
time limit, universal fluctuations of order is at a trans- 
verse scale y ~ ts, described by the scaling relation (1). 
We have probed numerically a refinement of this scaling 
law at short times, for the case of a random environment 
presenting correlations at a typical finite lengthscale £ 
and at finite temperature, by testing the scaling rela- 
tion (46) for the two-point free-energy correlation func- 
tion C^(t, y) defined in (16). If this scaling extends to the 
higher order correlation functions, it will be equivalent to 
stating that F(t, y) scales in distribution as 

i 

F(t, y) ~ [5 * X_s_ (vlVW)) ■ (49) 

Compared to the large time asymptotics (1), this rela- 
tion expresses the fact that to the transverse fluctuations 
(y 2 (t)) — B(t) correspond free-energy fluctuations of or- 
der J B(t), described by the roughness B(t) not only in 
the RM regime where B(t) oc fs but also at smaller times. 
The properly rescaled function Xg(y) would thus be in- 
dependent of the parameters c, D, T and £, the constant 
D capturing all the parameters dependence. We have 
in particular recovered numerically that the value [56] 
D = ^Sr describes correctly the high-temperature regime 
where the disorder correlation length £ plays no role and 
correlations take the form (10). Moreover, the distribu- 
tion of the variable XfjS) is known at £ = (xo(y) — x(y) 
being then the Airy 2 process [45]) and an interesting open 
question is thus to characterize the analogous process at 
non-zero £. One can conjecture that a suitable character- 
ization is provided by an appropriate Macdonald process 
(or a generalization of it) - already known [65, 79] to yield 
the process x(v) at £ = in some specific limit - espe- 
cially in view of recent results [80] allowing to represent 
Macdonald processes in terms of Brownian motions inter- 
acting within a finite range. This could in particular help 
to characterize the universality of the process Xi{y)i *- e - 
to determine how many details of the correlator (y) of 
the microscopic disorder (13) are left in x^(y). 

A point of particular interest is the scaling of fluctua- 
tions with respect to temperature. Regarding the rough- 
ness B(t) it has been shown [23, 52, 76] that, above a 
characteristic temperature T c = (£cD)a, B(t) can be 
rescaled as 



B(t;c,D,T)=^ h B(t/U 



with £ 



th 



cb 



1 



(50) 
(51) 



while below T c such a rescaling does not hold. Never- 
theless, our results on the free-energy correlator C(t,y) 
indicate that the extended scaling law (46) holds below 
and above T c at large enough times, the dependence in 



12 



the parameters (c, D, T, £) being then gathered in a sin- 
gle prefactor D. If the roughness exponent £rm is not 
expected to change below T c , however the prefactor A of 
the roughness in B{t) ~ Ai 2< > RM can be modified by £. 
This question is relevant in particular to determine the 
precise dependence of the crossover timescales of one- 
dimensional interfaces described by the DP coordinate 
y(t) [23, 53]. We refer the reader to [58] for numerical 
and analytical results in that direction. 

On the other hand, recent experiments on liquid crys- 
tals [42-44] (see also [9] for a burning front experi- 
ment) have demonstrated that the height fluctuations of 
a growth interface are correctly described by the (cen- 
tered) free-energy F(t, y) itself. In particular it was found 
that different Airy processes of the KPZ universality class 
successfully account for the observed height distribution 
(Airy! or Airy 2 depending on the geometry of the initial 
condition). The authors of those studies were in par- 
ticular able to measure the correlator C(t, y) with high 
precision. Given the existence of finite disorder correla- 
tion length £ in experiments, it would be interesting to 
study the equivalent of the low- and high-temperature 
regimes when tuning the system parameters which plays 
the role of T in our description. In the notation of [44] , 
the relation between the amplitude Y of x(y) (linked to 
our D in (49) by T = D 2 /c) and A and v in the KPZ 
equation (linked to our parameters by c = i- and T = K- ) 
would change from low v to high v. Precisely, defining the 
characteristic value v c — (£A 2 D) 3 , we predict a crossover 
from the (observable) dependence T ~ D 2 \v~ 2 at high 
v (y ^> v c ) to a saturated regime T ~ D 2 \v~ 2 at low v 
(y -c v c ). The influence of £ could thus be probed by 
measuring the power-law dependence of Y in the param- 
eter A displaying a non-trivial dependence Y ~ A~s for 
v <C v c . 
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Appendix A: Different formulations of the 
free-energy correlator 

We determine in this appendix some properties of the 
free-energy correlator C^(t,y) defined in (16). Defined 



by the decomposition (15) of the free-energy, the reduced 
free-energy F(t, y) is statistically invariant by translation 
along y, which implies that all its moments are invariant 
by translation along y at fixed time: 

F v (t, yi+Y)... F v {t, y n + Y)=F v (t, yi )... F v (t, y n ) . 

(Al) 

In particular, F(t, y) is independent of y, so that C^(t, y) 
is also the connected correlator of Fy(t, y). Indeed, start- 
ing from (15) 



v 2 y 2nTt 

Fv{t,y) = c y - + -\og +F v (t,y) (A2) 

At I C 



y 

c h cte(t) 

2t v ' 



thus 



y 2 T 2nTt 



F v (t,y) = F v (t,y) - --log 



F v (t,y)-F v (t,y) + cte(t) 



(A3) 



(A4) 
(A5) 



and hence from the definition (16), and dropping the in- 
dex V for simplicity 



C. 



:(t, y'-y) = { [F(t, y>) - F(t, y>)] - [F(t, y) - F%yj] }' 

(A6) 



[F(t, y>) - F{t, y)} 2 - [f(*, y') - F(t, y)] ' 



(A7) 



as announced. The fact that this expression depends 
only on y' — y is not obvious and arises from the STS 
through (15). See Ref. [58] for a study of the time de- 
pendence of the quantity Fy(t,y) which is independent 
of y. 

Note that the second derivative of C^(t, y) with respect 
to y is directly related to the correlator of the "phase" 
r] V (t,y) = d y F v {t,y). Indeed, using F v (t, y)-F v (t, 0) = 
Jo dyirjv(t,yi) one has 



d 2 y C^t,y) = dl[F{t,y)-F{t,Q) 
rv ry 



(A8) 



d y / dyi / dy 2 vv(t,yi)w(t,y 2 ) (A9) 
Jo Jo 
rv rv 



d 2 / d Vl / dy 2 R(t,y 2 -yi) (A10) 



= R(t,y)+R(t,-y) 
d 2 C ( (t,y)=2R(t,y), 

where the correlator of rjv is denoted by 



R{t,V% - 2/1 ) = ri v {t,y x )riv(t,y 2 ) 



(All) 
(A12) 



(A13) 
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Appendix B: Large y behavior of the 
correlator C$ (t, y) 

In this appendix, we show that, as discussed in sec- 
tion III A, the correlator C(t,y) presents 'wings' at all fi- 
nite times, i.e. goes to a finite constant limiyi-Hx, C(t,y) 
at large \y\ for t < oo. 

To this aim, let us first show that the inte- 
gral Af(t) — J dy R(t, y) of the two-point correlator of 
r]v(t,y), defined in (A13), is a conserved quantity. 
From (7), the evolution equation for the reduced free- 
energy Fy(t,y) of the directed polymer, defined in (15), 
is given by (43) 

dtFv(t,y) = ^-fiFv^v) ~ ^[dyF v (t,v)] 2 

~^d v F v (t,y) + V(t,y) (Bl) 



so that the evolution equation of rjv(t,y) = d y Fy(t,y) 
also includes an explicit time dependence: 



dtvv{t,y) =—^d%r]v(t,y) ~ -r)v(t,y)d y ri v {t,y) 

1 y 

- TVv(t,y) - -dyVv{t,y) + d y V(t,y) . 

(B2) 



One can now compute dt J dyR(t,y). Taking advantage 
of the invariance of the disorder distribution through the 
symmetry by reflection y i— > — y, and of the statistical 
invariance of r]v(t,y) by translation along y, one has, 
dropping the index V for simplicity 



dt / dyR{t,y) = d t / dyy(t,y) V (t,0) 



I 



dyly(t, y) 0) - ~t?(t, 0)d y y(t, 0) - l v (t, 0) + 8 y V(t, 0) 



d t / dyR{t,y) = -- / dyR(t,y) 



(B3) 
(B4) 



+ : 



V(t, 0) [jfrftTlfo V) - -»?(*. y)dyV(t, y) —-*?(«, y) - ^d yV (t, y) +8 y V(t, y)] 



-jd y [yr/(t,y)} 



dy T)(t, y) V (t, 0) + —ri(t, 0) [d yV (t, y) - d yV (t, 0)] y y ZiZ 



^(*,o)M*,v) a ]^:-^(* I o)[ W (*,v)]^ 



(B5) 
(B6) 



This shows that Af(i) verifies the equation dtM(t) = 
— jAf(t) whose solution reads Af(t) — % . The constant 
C\ is equal to thanks to the initial condition Af(0) = 
(the initial condition for F is -FV(0, y) — 0, see part III A, 
and thus rjv{0, y) = 0). This yields that Af(t) is constant 
in time and equal to 0: 

Af(t) = J dyR(t,y) = 0. (B7) 

This equality at all finite times t is equivalent to the 
existence of wings in the correlator C^(t, y) at finite times. 
Indeed, since 2 / Q y dy' R(t,y') = d y C ( {t,y) [see (A12)], 
this equality also writes 



lim d y C 6 (t,y) = 0, 



and this corresponds to the plateau at large \y \ of Fig. 3. 
This result is already known to hold at large finite time 
and for a delta-correlated disorder (£ = 0) [45]. In this 
appendix we thus have generalized this result at all times, 
and to a disorder correlator with short-range correlations. 



Appendix C: Short-time dynamics (diffusive scaling) 

In this appendix, we study in the short-time regime 
and at finite £ the scaling behavior of the two-point cor- 
relators C^(t,y) and R^(t,y) of Fv(t,y) and r)v{t,y) — 
d v Fv{t,y), defined respectively in (16) and (A13). 

The polymer is pinned in y = at time t — so that 
the initial condition translates for the reduced free-energy 
into Fy(0,y) — 0. Let us thus assume that Fv(0,y) 
remains small at short times so that in this regime the 
non- linear term of (Bl) remains small and negligible com- 
pared to the linear terms: 



(B8) d t F v (t,y) 



Y c 9yFv(t,v)- 



d y F v (t,y) + V(t,y). (CI) 



This equation is linear and can thus be solved directly. 
Before doing so, note that for an uncorrelated disorder 
(£ = 0), the corresponding steady-state distribution of F 
is the same [56] as in the non-linearized one (9) (remark 
from (15) that F and F share the same distribution at 
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infinite time). In particular, denoting C^ n (t,y) (resp. 
R^ n (t,y)) the same correlator as (16) (resp. (A13)) but 
for F v solution of the linearized equation (CI), one also 
has 



C l t(t,y)k= 



cD 
-Tfrlvl 

> — %) 



it— >oo ± 



(C2) 
(C3) 



where i^ in (t,j/) = \d>*Cf>-{t,y) as in (A12). Although 
the solution of the linearized evolution equation (CI) is 
valid only at short time compared to the complete evo- 
lution (Bl), it is thus instructive to study its behavior at 
all times. 

In what follows, we determine the finite t and finite £ 
equivalents of the correlators (C2,C3). To get rid of the 
term \d y Fy in (CI) we set 



F v (t,y) = Vie^F v (t,y) 
The evolution of Fy is then 

d t F v (t,y) 



(C4) 



V(t,y) 



Besides, in absence of disorder (V = 0), 



G(t,y) 



2nt 



6{t) 



(C6) 



is a Green function of the equation for Fy, that is: 

[0t - ^dl] G(t -t',y- y') = S(t - t')S(y - y') . (C7) 

We denote by 6(t) the Heaviside step function. The so- 
lution of (C5) is hence 



dt' dy'G{t-t\y-y')V(t\y'). (C8) 



F(t,y) 



Here and in what follows, integrals along direction y run 
by convention on the real line. 

Let us first determine the effect of £ by computing the 
value of the peak of the correlator 



1 



Rt(t,0) = ^[F(t,y)-F(t,y)y 



v=o 



(C9) 



= d y F(t, 0) d y F(t, 0) 



which is a matter of Gaussian integration. One has 
d y F(t, 0) = Vtd y F(t, 0) and 



+00 



d v F(t,0)= I dt 1 \dy l d y G{t-t 1 ,-y 1 )V{tx,yi) (C10) 
I /", /', ,■ f3c \1 ggj 



V2ir Jo 



dt 



<///:(, ', . ) !l " J " 7 



-V(h, yi ) (Cll) 



so that putting everything together, and considering for 
the disorder V the Gaussian corr elator Rf mss (y) defined 
in (14), one gets 



i?£ n (t, 0) = t d y F{t, Q)d v F(t, 0) 



t ldy 1 dy 2 I dt 1 dt 2 ^ " " 



2tt 

(/?c) 3 



£c / ^1 | ^2 
2 V t — ti ~ l_ t — t 



i *-*2 / e 



0c f«? 2 ^ 



(1/2— Vl) 

4 «" 5(ti-t 2 ) 
. 



(<-ti)i(t-£ 2 )2 



V(yi,ti)V(y 2 ,t 2 



4tt5^ 7 Jo (t-hYh 



ViV2 



exp 



2 vt-" 7 ; + ^J (yi+y2) ~^^ 



/3c 



(C12) 
(C13) 

(C14) 



The quadratic form in the last exponential writes and finally (C14) becomes 



— ^y T Ay with y = (2/1,2/2) and A the matrix 



A = 



1 



1 



(C15) 

The result of the Gaussian integral with respect to y is 



2tt 
^4 



{/3ct^ 2 t + ii(t-fi)] } 



3/2 



(C16) 



: PCD 



3/2 



4tt2^ Jo {^ct[^ 2 t + ti(t-ti)]} 



t(t + 2/?c£ 2 ) - £y/Pci{t + 4 / Se£ 2 ) arccot (2£ J^) 



(C17) 



2^C(i + 4/3cC 2 



(C18) 



We see on one hand that the infinite time limit reads 

cD 



limi? "(<,()) 



(C19) 
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which diverges indeed as £ — > 0. The effect of £ is to regu- 
larize the £ = result (C3) around y = 0. The factor 
directly arises from the Gaussian correlator (14) as will 
become clear below. On the other hand, the short-time 
behavior writes 



R [m (t,0) 



Dt 



12£ 2 



0(i 3 / 2 ) 



(C20) 



As for the correlator C^ n (t, y) of the reduced free-energy, 
this yields from (A12) the small y expansions at large 
and short times 



that the reduced free-energy Fy (t, y) of the polymer end- 
point (the 'effective potential' it sees) has a Gaussian 
distribution, i.e. that it is fully determined by its mean 
value and its two-point correlator C toy (t,y). 

Following the scheme used in Ref. [53], we com- 
pute in this appendix the Gaussian- Variational-Method 
(GVM) approximation of a DP toy model with satura- 
tion 'wings' appearing at \y\ > It at time t. We check 
specifically that by consistency they should appear at 
it ~ y/B(t) at asymptotically large lengthscales, as as- 
serted in Sec. IIIC. 



C l t(t,y) 



C x t(t,y) 



cD 



(t->0) Dt 



y for y <C 



for y <C 



(C21) 
(C22) 



The full (t, y)-scaling of the correlator C^ n (t, y) can be 
determined using the same approach as exposed above, 
starting from the solution (C8) for Fy(t,y). One finds 



q*(t,y) = D ( & x - ^~ 
Jo 



'7T 



-e 



' ti(t-ti) 

/3ct 

(C23) 

To reveal the scaling of this expression, one performs the 
change of variable ii = tr in the integral 



cD 



1 1 



eXP L 4 ({^ C /t + T(l-T))J 



/3c «V^7* 



(l,y^Mt) 



(C24) 



which takes the scaling form, with the thermal roughness 
-Bth(i) = ^ = " 



where the scaling function is 

Cf(y) = C f (i = l,y)| c=/3=D=1 
exp [ 



i 1 
dr- 



4(« 2 +r(l-r))J 



yfty/T(l - T) + £ 2 



(C25) 

(C26) 
(C27) 



The relation (C25) describes a scaling form of the free- 
energy correlation, with the roughness B t h {t) of the ther- 
mal regime. We discuss in sections IV B and VB an ex- 
tension of its validity to the large time regime of the DP, 
with the full roughness B(t) instead of B t -^(t). 



Appendix D: Solution of the winged DP toy model 
in the Gaussian Variational Method approximation 

As introduced in Ref. [71, 81], what we call generically 
a DP 'toy model' is essentially based on the assumption 



1. Replicae 

The replica approach allows to determine the statisti- 
cal average of an observable 0(t, y) of the DP toy model 
as defined in paragraph III C from the formal expression 



(0(t,y))= Km / d yi ...dy n O( yi ) e -P£S=i *■"*<*.».) 

(Dl) 

= lim / dyi...dy n O(y 1 )exp{-0F(t,y)} 

(D2) 

where the replicated free-energy F(t,y) of n copies 
y = (j/i, . . . , y n ) of the polymer endpoint reads 



Tl Q 71 

" + 7 E C to y(t,y b -y a ) . (D3) 




„b=l 



F di .(t,y) 



To obtain this expression, we have used the assumed 
Gaussianity of the F toy distribution, presumably inher- 
ited from the random-potential Gaussian distribution. 



2. Gaussian Variational Method (GVM) 

There is no known way of determining the exact rough- 
ness of the toy model. The GVM approach [71, 81] for 
the toy model (see [53] for a detailed presentation) con- 
sists in approximating the replicated free-energy (D3) by 
a trial quadratic free-energy 



a,b=l 



(D4) 



parametrized by the n x n hierarchical matrix G 1 (t): 
G- b 1 (t) = C -5 ab -a ab (t) (D5) 

with the connected part Gj 1 (i) = J2 b ^ab = f nxe d by 
the case in absence of disorder. The corresponding rough- 
ness is directly read from the diagonal term of its inverse 
matrix, combining (D2) and (D4): 

B(t) = T lim Gaa(t) . (D6) 

n— !-0 



16 



To find the best quadratic approximation of F(t,y), 
the extremalization conditions read for the pairs {a, b}: 



dJ v; 



dG ab {t) 



[G(t) G -1 (t)] =0 



(D7) 



where the variational physical free-energy J var (i.e. after 
averaging over thermal fluctuations) is defined at each 
time t as 

Jvar = T Q + (F - F )o = Jo + < F el - F )o + (^dis)o (D8) 
and the trial physical free-energy is 

J = —T log Z = - - log dct G + Cte (D9) 

In those expression, (-)o denotes the average with respect 
to the normalized Boltzmann weight e~ /3F ° /Z and Z is 
the corresponding partition function. Let us first com- 
pute 

n n 

(for4 = (^^,E VaG-J{t)y b ) (D10) 



a=l 



lj2[ C t S ^- G ab(t)](Vay b )0 (Dll) 
ab 



cT 
2i 



Y,Gaa(t)+ Cte 



(D13) 



and using the Fourier transform of the correlator 
C* toy (i,y) as defined in (25): 



(F dia ) 



E 



dX 2 
2^ A 5 



<l-cos [\(y b -y a )]) B*°*(t,X) 



(D14) 



pB^ [dx i- e -^ 2 <^^) 2 >o 



2 



E 



2^ Y 2 



(D15) 
Rf(t,X). 
(D16) 



Gathering the results (D9)-(D13)-(D16) in the defini- 
tion of J va r(£) (D8), we obtain an explicit expression of 
Jvar [G(i)] and we can apply the extremalization condi- 
tion (D7) for the off-diagonal terms a ^ b (setting the 
usual notation G = G aa )- 

(D17) 
(D18) 



since the sum ^ ab in (D16) contains twice G a b by sym- 
metry of the matrix, hence the '2x' in (D17). Thus 
for a continuous parametrization a(u) (u £ [0,1], see 
Ref. [82]) of the hierarchical matrix a a b, in a full replica- 
symmetry-breaking (full-RSB) formulation, one gets the 
self-consistent equation at each time t: 



a(t,u)=2t 



D f dX 



-\ 2 T[G(t)-G{t,u)] 



R~ oy (t,X) (D19) 



where [G(t),G(t,uj\ characterizes the matrix G(t) 
whose matrix inverse has off-diagonal coefficients 
G- 1 (t,y) = -cr(t,u). 



3. Formal solution of the equation for a(u) 

The algebra of hierarchical matrices in the limit n — > 
has been worked out in Ref. [82] , and following the con- 
ventions of Ref. [53] we recall the properties needed there- 
after in the derivation of the GVM solution: 



d u [G-G(u)] =- 



[a] (u) 



UG 



G - G(u) 



(g„ 1 + [ct](«) 

(u) - J^dva{v] 
f 1 i. 

Ju V 



1 f 1 dv 

2 r^-\ 



u G- L + [a]{u) 

1 dv [cr}(v) 



G-' + [a](v) 
G- 1 



(D20) 

(D21) 
(D22) 

(D23) 



r - -J— 1 1 4- r <k 

~ G7 1 \ + J « 2 Gc l + [<r](v 

where the connected term G~ 1 (t) = c/t. Denoting 

Jk(t)= ( — X k e~ x2T ^~ G ^Ri oy (t,X) (D24) 

the extremalization equation (D19) becomes simply 
<j(u) = 2$DJq. Differentiating with d u and assuming 
that <j'(u) 7^ 0, we obtain 



2D 



[G- c l + [a](u)] 



>J2 



(D25) 



Isolating G^ 1 + [cr](u) = (2DJ 2 ) 1 ^ 2 and differentiating 
again with d u yields: 



T\/ 2D 



[Go 1 + [a]{u)\ 



(D26) 



One finally obtains that a(t, u) can be determined by 
solving the coupled self-consistent equations 



a(u) = 2pDJ and u = 2^ 2 D' 2 T 



(D27) 



T 

= -^Pab - F> 



2tt « 



and the definition (D24) which hides the ^-dependence, 
e.g. by isolating G — G(u) as a function of u from the 
second equation, and then inferring a(u) from the first. 
These equations are valid on segments of u £ [0,1] where 
u(u) is not constant. 
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4. Solution for Cauchy wings at large times 



Up to now, the results are valid for a generic toy model 
correlator i?~ oy (i, A). In Ref. [53] we have worked out the 

case of a rounded correlator i?~ oy (t, A) = e~ A "^~ . In order 

to add saturation 'wings' to this correlator, we consider 
a Cauchy Ansatz (see paragraph III C and Fig. 3) : 



X 2 



-A £ 



(D28) 



where it corresponds to the transverse lengthscale at 
which the wings start to develop. From the infinite-time 
result (18), we expect i t to diverge as t goes to infin- 
ity, which allows us to study the large scale regime by 
focusing on a small i^ 1 expansion. However, bluntly ex- 
panding i?| oy (i, A) in powers of e^ 1 < 1 in (D19) yields 

divergent integrals; an alternative solution is rather to 
solve the coupled equations (D27) perturbatively in . 
Expanding Ji and JZ4 in powers of i^ 1 one finds that the 
following equation in X: 



<© 5 



(D29) 



admits X = £ 2 + T G(t) — G(t, u) as a solution to the 

first order in . The physical solution must go in the 
limit ^ to the solution obtained in Ref. [53] at 
l^ 1 = 0, which involves two RSB cutoffs it*(i) < U c (£) 
such that the solution a(t, u) is non-constant only for 
it*(i) < u < it c (£) and times below the Larkin time t c de- 
fined by u*(i c ) = u c (£). We have actually: 



G(t)-G(t,u) 



3 4 ?rT 4 . 3 7 7tt 2 T 8 „_ l , 



2 8 D 2 



2 W DH 



(D30) 

Inserting this result into the first relation in (D27), one 
finds that 



a(t, u) 



2D 2 flu 
k^F I 3T 



+ 0(£ t " 2 ) 



(D31) 



whenever a(t,u) is not constant. To minimal order in 

l^ 1 , a(t, u) is translated by the constant — ^£^~. We also 
see from its definition (D21) that [a](t, u) is not modified 
at this order: 



[a](t,u) = [a](t,u)\ e . 1=Q + 0(^ 2 ) 



(D32) 



Besides, the definition (D21) imposes the condition 
[<r](t, 0) = and thus the existence of u*(t) for which 



[a](t,u< u*(t)) = 



(D33) 



Both a(t, u) and [a] (t, u) are constant in u below (t) 
and can be strictly monotonous only above w*(t); this 
implies by continuity of a(t, u) that 



a{t,0) = <r(t,u< u»(t)) 



D 



(D34) 



To determine whether the second threshold u c (£) has a 
correction which depends on i^ 1 at this order, one has to 
check that the original equation (D19) is satisfied. Using 
the replica inversion formula (D22), one checks that u c 
being fixed this quantity is not modified at order 0(£^ 1 ): 



G(t)-G(t,u) = [G(t)-G(t,u)] 



■0{tl 2 ) (D35) 



and replacing its value in the original self-consistent equa- 
tion for a(t,u) (D19) we see that the equation on zt c (£) 
is not modified at the same order. 

We can hence compute the modification to the rough- 
ness B{t) = T\im n - j .QG(t) as given in (D6) through the 
formula (D23). Observing that 

<r(t,0)-a(t,Q)\,-i n 

°">' c <"i..-- + [C A,]- +o(c) 

(D36) 
(D37) 

(D38) 

(D39) 



G(t)\^ =a - t - 2 ^7 1 + 0(ii 2 ) 



we finally obtain 
B(t)=B(t)\ 



t 2 D 



l =o 



^r 1 



At asymptotically large times, we know[53] that 
_ 3 f2D 2 - 



B(t) 



17 =0 



3,1 
t3 



2 \ TTC 4 

so in the random-manifold roughness regime 

_ 1 

BR M {t) ~ ff^rj 3 t^- Thus for the lengthscale £ t 
not to modify the scaling of this asymptotic behavior 
via the correction — t —§-l~^ x 1 we must have consistently 



Appendix E: From the discrete to the continuous DP 

In this appendix, we determine formally the contin- 
uum limit of a discrete directed polymer model, whose 
parameters are denoted for convenience in Gothic script 
(the inverse temperature ft, the disorder strength 2), and 
when needed the elastic constant c). See also Ref. [52] 
for a detailed analysis of the scalings and Ref. [61] for a 
mathematical approach. 
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1. Simplest case: jump to the two nearest 
neighbors (n = 1) 

The partition function Z ty of the discrete SOS model, 
where the polymer can jump one step either to its right or 
to its left at each time step, obeys the recursion relation 



+ 1 



(El) 



It is described by an inverse temperature ft and a disorder 
strength 35 (hidden in V through Vt,yVt' ,y' = QStt'Syy 1 )- 
To explicit the correspondence between the discrete pa- 
rameters ft and 25 and the continuum parameters /3, c 
and D (e.g. of the evolution equation (6)), we may ex- 
plicit the lattice spacings a and b in directions y and t 
respectively: 



Zt-b,y-a + Zt-b.y+a 



(E2) 



where the factor (abD) ? in front of the disorder is chosen 
so that V^ y = (abD)~^Vt, y becomes a white noise in the 
continuum limit: 



1 

abD 

a,b-¥0 



b a y 



>5(t-t')S(y-y') 



(E3) 
(E4) 



Introducing for normalization purposes 
W t .y = 2~ t Z t , y , one expands its corresponding equation 
of evolution as follows: 



W t -b,y-a + W t - b ,y+a 



(E5) 



~[l-ft{ab®)Wt v ] \2W t - blV + a 2 d 2 y W t J (E6) 



so that, with W t>y — W t -b,y — bd t W t . y 



9tW t , y = ~tfW t , v - ft{ab- x ®)ivl y W t , y (E7) d t W t , y = 



2 b 



which corresponds for instance to the continuum 
Feynman-Kac evolution (6) with the parameters 



(3=fl c = 



b 1 

a? ft 



D = a6 _1 D (E8) 



This correspondence between discrete and continuum pa- 
rameters is used in section IV when discussing numerical 
results on the discrete DP. 



2. Generic case: jump to the 2n nearest neighbors 

For the sake of completeness, we now consider a gen- 
eralized SOS model in the spirit of Ref. [52], where the 
polymer endpoint y can jump to either of its 2n neigh- 
bors y — 2n+ 1, y — 2n + 3, . . . , y + 2n — 1 with an elastic 
weight depending on the distance j as e - -^ . The parti- 
tion function now obeys 



Zt,-i 



3=1 



Zt-l,y-(2j-l) + Zt-l, y +{2j-l) 



(E9) 

We explicit as previously the lattice spacings a and b, 
and we also introduce a microscopic elastic constant c 



Z t ^ y = e-^ ab ^v t ) yx 



3=1 



(E10) 



J t-b, y -(2j— l)o + Zt-b,y+(2j-l)a 



Introducing, for normalization purposes, Wt, y — 
(2n)~ t Z t y one expands its corresponding equation of 
evolution as follows: 



W t ,y~^[l-ft(ab<3)^Vl y ]x 



(Ell) 



E e ~ 4, 

i=i 



2W t -b, y + {2j - lY^d'Wt-b 



To ensure that the dominant order of the right hand side 
is Wt—b,y, which allows to recognize a time difference, one 
sets 



" = E e ~ 

3=1 



(E12) 



Then, defining the effective elastic constant n n (c,ft) as 
in Ref. [52] 



En 
.1=1 



e 2 



£Ei=i(2j-l) 2 e-5 co2 « 2 



(E13) 



one recovers from (Ell) the equation of evolution 
1 a 2 



2p Kn b"y v,t ' v 



%W t ,y - fi(ab- 1 l))tWt,yV f , v (E14) 



which corresponds to the continuum evolution (6) upon 
the identification 



P = fi 



;K„(c,fi) D = ab~ x, D (E15) 



Note that the result holds for any n; in particular, K n = j 
for n=l (as obtained in the previous paragraph) but K n 
goes to another c-dependent limit for n 2 3> T, which 
enables in particular to study the high-temperature limit 
of the discrete directed polymer model [52]. 
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